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^-H , Abstract. A multivariate Gauss-Lucas theorem is proved, sharpening and gen- 

erahzing previous resuhs on this topic. The theorem is stated in terms of a 



a^ 



^ 



seemingly new notion of convexity. Applications to multivariate stable polyno- 



^ I mials are given. 
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t^ . Let P{z) be a polynomial in the complex variable z with non-null derivative Q{z). The 

classical Gauss-Lucas theorem is the relation 

^ . g-^(O) C H{p-\0)) (1) 

CN I where H{P ^(0)) is the convex hull of the roots of P{z) in the complex plane. (See Marden 

^ . [5, Theorem 6.1].) It is intriguing to look for a multidimensional extension of this elegant 

^ ■ result. 

O ■ 

CN ' To establish notation, let P{z) stand for a polynomial in the M complex variables z = 

{zi, . . . , zm)- Given P{z), let Q{z) stand for the partial derivative of P{z) with respect to 

Zk- The desired extension of ([1]) has form 



Ql\Q) C H{P-\Q)) , (2) 
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under the condition that Qk{z) is non-null. The convex hull in (|2]) is taken with respect to 
the usual convex structure on C^^ treated as isomorphic to i?^*^. (C stands for the complex 
numbers and R stands for the real numbers.) 

It is somewhat curious that no proof of ([2]) seems to be available in the literature under 
the conditions stated. The papers [3] and [4] prove a direct generalization of ([1]) in the form 

M 

fl Q-,\0) C H{p-\0)) , (3) 

k=l 

under the condition that P{z) has finitely many zeros. (The left-hand side of ([3]) represents 
the critical points of P{z) in C^^.) 

The finiteness restriction on the null set oi P(z) in ([3]) is removed in a related result 
about multivariate stable polynomials, i.e. polynomials that have no roots z= {zi, . . . , zm) 
with Im{zk) > for all k. It is shown in [1, Theorem 3.1] and [7, Lemma 2.4] that the 
k-th partial derivative Qk{z) of such a polynomial P{z) is also stable, a result called the 
"multivariate Gauss-Lucas theorem", even though the roots of Qk{z) are not shown to lie 
in H{P~^{0)). This further refinement is provided in Theorem 1, using a weaker notion of 
convexity than that used to define //(P~^(0)), and thereby obtaining a smaller convex hull 
for P^\0). 

Preparatory to the proof of Theorem 1, it is necessary to establish some terminology. For 
z G C*^ with z={zi, . . . , Zm), the variables Zk are called the coordinates of z. The notation 
z^'^'^ stands for the M — 1 coordinates of z obtained by omitting the coordinate Zk- Given a 
subset A of C*^, the section of A determined by z^''^ is defined by 

A{z')^{ykeC:yeA, y^''> = z^''>} . 

A set A C C^^ is called separately convex in C^^ if A{z^''^) is a convex subset of the 
complex plane for all k G {!,... ,M} and all z G C^^. If two sets are separately convex 
in C*^, then so is their intersection. This allows the definition of the convex hull of a set 
A C C*^ as the smallest set which is separately convex in C^^ and contains A, denoted 
H2iA). 

Theorem 1 If P{z) is a multivariate polynomial defined on C^^ , with non-null k-th partial 
derivative Qk{z), then 

Q,\0) C H,{p-\0)) . (4) 



Proof. Without loss of generality let k = l. Suppose z G C^ with Qi{z) = 0. To prove 
the theorem it suffices to show that z G H2{P^^{0)). 

Let f{w) = P{w,z^^^) for w E C. Regard ,f{w) as a polynomial in w whose coefficients 
depend on the fixed vector z^^\ Note that 

by assumption. The univariate Gauss-Lucas theorem implies that zi G /J(/~^(0)), where 
H{f^^{0)) stands for the usual convex hull in C of the roots of f{w). Let A be any subset 
of C^ which is separately convex in C^ and contains P~^{0). Then the section A{z^-^'>) is 
a convex subset of C containing /~^(0); this implies that A(z^^') contains H(f^^(0)). In 
particular, Zi G A(z^^^), hence z G A. By choice of A, this shows that z G i^2(-P HO)) ^.nd 
proves the theorem. D 

The application of Theorem 1 to multivariate stable polynomials can easily be broadened 
by generalizing the definition of stability. Let 9 = {9i, . . . , 9m) be a real vector and define 

A{9) = {zeC^ : Im{e'^' Zk) > ior 1 < k < M} . 

Let A{9) stand for the complement of A{9) in C'^ . 

Lemma 1 A{9) is separately convex in C^^ . 

Proof. It suffices to prove that for z G C^ the section of A{9) determined by z^'^^ is a 
convex subset of the complex plane for all k. Without loss of generality set k = l and let 
Zj = Xj + iyj, where Xj and yj are real. Let c stand for the minimum of Im{e^^^ Zj) for j > 2. 
The section of A{9) determined by z^-^'^ is the set 

{x + iy E C : {x sin 9i + y cos ^i) A c < 0} . 

It is easy to see that this set is convex by considering the cases c > and c < separately. 

D 

Theorem 2 Suppose the multivariate polynomial P{z) is 9-stable, meaning that it has no 
zero in A{9). Then any non-null partial derivative Qk{z) is also 9-stable. 



Proof. Apply Lemma 1 to get H2{P ^(0)) C ^4(6') and finish using Theorem 1. D 

It is clear that i/al^^HO)) C //(P"^(0)), hence (HD sharpens ([2]). With a view towards 
further sharpening ([2]), it is natural to compare separate convexity in C^'^ with the analogous 
notion of separate convexity in R^ . The latter notion can be defined by substituting R for 
C in the definition of separate convexity in C^ and has been studied under many names. 
(See [2] and [6].) 

If C^^ is identified with R^^^ via the real and imaginary parts of the complex coordinates 
in C^^, then any separately convex subset of C^^ is separately convex in R^^ . Letting 
Hi{P~^{Q)) stand for the smallest separately convex subset of R^^ containing the roots of 
P{z), it follows that 

H,{p-'m c H,{p~'m . (5) 

However, (jl) is not valid if ifi(P"^(0)) is substituted for i72(-P~^(0)), even if M = 1 (Example 
1). Thus, no straightforward sharpening of (jl]) seems to be available. 

Example 1. When M = 1, the univariate case, conditions are given for cubic and 
quadratic polynomials P{z) so that (jl]) is not valid iiH2{P^^{Q)) is replaced with i^i(P^^(0)). 

If P{z) is a cubic polynomial with real coefficients and real roots, then ifi(P^^(0)) is 
equal to i^2(-P^^(0)), so the substitution does not affect the validity of (jlj). Therefore, assume 
that P{z) has roots a+hi, a— hi, and c (with a, 6, c real and & 7^ 0). Setting the coefficient of 
z^ in P{z) to be 1, yields that 

P{z) = {z-c){{z-af + h'') . 

The derivative Q{z) = P'{z) is calculated to be 

Q{z) = 3z^ - (4a + 2c)z + a^ + b^ + 2ac . 

The set iJi(P^^(0)) is connected, consisting of the line segment joining a + bi to a — bi and 
the line segment joining a to c. 

The roots of Q{z) can be written as 

w, = f + l±'-{3b'~{a-cr)K 

where 3b^ > (a — c)^ by the assumption that not all roots are real. If a 7^ c, then it is 
clear that the roots of Q{z) do not lie in the set //i(P^^(0)). If a = c, then i^i(P~^(0)) 



consists of the line segment joining a + bi to a — bi and the roots of Q{z) are of the form 
a±{^)beH,iP-\0)). 

In summary, if P{z) is a cubic polynomial with real coefficients, then the roots of Q{z) 
lie in Hi{P^^{0)) if and only if the roots of P{z) lie in a straight line parallel to either the 
real or imaginary axis in the complex plane. If P{z) is a quadratic polynomial with complex 
coefficients, then the same result is true. (Note that in the quadratic case, Hi{P^^{0)) 
consists of the two roots of P{z) if these roots do not lie on a line parallel to either the real 
or imaginary axis. Thus i/i(P~^(0)) is not connected in this case.) 
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